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Abstract. 

The total electronic Raman scattering spectrum, including the nonresonant, mixed and resonant components, is determined 
for the charge-density-wave (CDW) phase of the spinless Falicov-Kimball model at half filling within dynamical mean-field 
theory. Its frequency dependence is investigated for different values of the energy of the incident photons. The spectra reflect 
the different structures in the density of states and how they are modified by screening and resonance effects. The calculations 
are performed for the B\„, i?2g and Ai g symmetries (which are typically examined in experiment). Our results for the resonance 
effects of the Raman spectra, found by tuning the energy of the incident photons, give information about the many-body charge 
dynamics of the CDW-ordered phase. 
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INTRODUCTION 

Experiments on inelastic light scattering are employed to learn about the complicated charge dynamics of a wide 
class of strongly correlated electronic materials. The photon couples to the charge excitations during the inelastic 
scattering process and directly probes the charge excitations of different symmetries. In this work, we study the 
strongly correlated electron systems with charge-density-wave (CDW) ordering. CDW systems possess a static spatial 
modulation of the electronic charge with some spatial ordering wavevector. Since the underlying ionic cores are 
charged, they also respond to this charge modulation of the electron density and often create a distorted lattice structure 
that follows the modulated charge order of the electrons. A direct measurement of the lattice distortion due to the ionic 
displacement is often the best way to measure the presence of CDW order; it is more difficult to directly measure the 
electronic charge modulation in the material. 

In the present work, we investigate how CDW order affects inelastic light scattering experiments by utilizing 
dynamical mean-field theory (DMFT) to exactly solve for the total electronic Raman spectra. Since inelastic Raman 
scattering is sensitive to different symmetry charge modulations (when polarizers are used on the incident and scattered 
light) it can provide information about the symmetry of the CDW state and of the many-body charge excitations. We 
expect our results should be relevant to different experimental systems that display charge-density-wave order via 
nesting on a bipartite lattice at half filling, especially in compounds which are three-dimensional like BaBiC<3 and 
Bai_ v K t Bi03 [Q], B [3|] because DMFT is most accurate in higher spatial dimensions; it may also be relevant to some 
layered two-dimensional systems, at least in a semi-quantitative fashion. Our work is the next step in recent results 
on transport, optical conductivity, and nonresonant inelastic X-ray scattering in CDW systems Jj,|5l|a] to the realm 
of resonant inelastic light scattering. Since experimental inelastic light scattering work on CDW systems was focused 
on the Raman scattering of the soft phonon modes, future experimental work should examine the electronic scattering 
directly. Hence this work has the potential to be directly relevant to the next generation of Raman experiments on 
strongly correlated CDW materials. 

One of the simplest models which posseses static CDW ordering at low temperature is the Falicov-Kimball 
model [7]. Historically, this model was introduced in 1969 to describe metal-insulator transitions in rare-earth com- 
pounds and transition-metal oxides. Later, it was found that it has an exact solution within DMFT H (for a review see 
Ref. [9]). The Falicov-Kimball model has two kinds of particles: itinerant electrons and localized electrons. Mobile 
electrons hop from site to site with a hopping integral between nearest neighbors and they interact with the localized 
electrons when both sit on the same site (the interaction energy is U); we denote the itinerant electron creation (an- 



nihilation) operator at site i by dj (</ ; ) and the local electron creation (annihilation) operator at site i by f t (f t ). The 
model has commensurate (chessboard) CDW order at half filling and this is the main property we exploit here. Brandt 
and Mielsch determined the formalism for calculating the ordered-phase Green's functions II 1 QTI shortly after Metzner 
and Vollhardt introduced the idea of the many-body problem simpli fying in large dimensions fi"ill . The CDW order 
parameter is known to display anomalous behavior at weak coupling lUHt], and higher-period ordered phases exist 
on the Bethe lattice Q. 



DMFT FOR THE CDW PHASE OF THE FALICOV-KIMBALL MODEL 

As mentioned above, the Falicov-Kimball model possesses the possibility for a transition into a commensurate CDW 
phase with doubly modulated (chessboard-like) density of charge, when both the itinerant and localized particles are 
half-filled. Since the hypercubic lattice is a bipartite lattice, implying that it is separated into two sublattices (called 
A and B) with the nearest-neighbor hopping being nonzero only between the different sublattices, the CDW order 
corresponds to the case where the average filling of the electrons remains uniform on each sublattice, but changes 
from one sublattice to another. We start by writing the Falicov-Kimball model Hamiltonian as the sum of its local and 
nonlocal parts 

ia ijab 

where i and a = A or B are the site and sublattice indices, respectively, and is the hopping matrix, which is nonzero 
only between different sublattices (tfj = tfj B = 0). The local Hamiltonian is equal to 

Hi =Un id n if - \l d n id - \i f n if ; (2) 

with the number operators of the itinerant and localized electrons given by n,v/ = d\ d i and «y = fjfj, respectively. Note 
that we have introduced different chemical potentials for the different sublattices. This is convenient for computations, 
because it allows us to work with a fixed order parameter, rather than iterating the DMFT equations to determine the 
order parameter. That method of iterative solution is subject to critical slowing down near T c , while working with a 
fixed order parameter is not. The equilibrium solution occurs when the chemical potential is uniform throughout the 
system (p£ = fijj and = jU?), which is the unique condition used to find the order parameter at a given temperature. 

We apply the DMFT, which provides an exact solution for the Falicov-Kimball model in the limit of infinite spatial 
dimensions. In contrast to the uniform case fl5l[l6ll . in the CDW phase, the DMFT equations become matrix equations. 
Since the DMFT solutions for the chessboard phase are described in detail in previous work we concentrate 

only on a few basic points as a summary and to establish our notation. The first step of DMFT is to scale the hopping 
matrix element as t = t* /2a/D flUl (we use t* = 1 as the unit of energy) and then take the limit of infinite dimensions 
D — > oo. The self-energy is then local: 

= £ fl ((»)^S a6 , (3) 

and in the case of two sublattices has two values T. A (co) and T. B (co). Now, we can write the solution of the Dyson 
equation (in momentum space) in a matrix form 

G k ((o) = [z(co)-t k }- 1 , (4) 
where the irreducible part z(o) and hopping term are represented by the following 2x2 matrices: 

z(co) -( » + \ U -( Q £k \ (5) 

with e k = t * Hm Y$=\ cos(ki)/VD. 

The second step of DMFT is to map the lattice Green's function onto a local problem by means of the dynamical 
mean field. Since there are two sublattices, a dynamical mean field X"(co) is introduced on each of them. As a result, 
the local lattice Green's function on each sublattice becomes: 

Gaa(C0)= co + n»-L«(co)-X°(coy (6) 



The third equation that closes the system of equations for G aa (co), E"(o) and X a (co) is obtained from the condition 
that the local Green's function is defined as the Green's function of an impurity problem with the same dynamical 
mean field X"(co). Such a problem can be exactly solved for the Falicov-Kimball model and the result is equal to 

1 - n a f n a f 
C aa (m\ — - I - n\ 

where n c j is an average concentration of the localized electrons on sublattice a which is found from the equilibrium 

condition of a uniform chemical potential (p.^ — fl 1 ^ — 0). 

These equations are self-consistently solved numerically. In Ref. [5], we analyzed the evolution of the DOS in 
the CDW-ordered phase. We summarize the main points which are needed here. At T — 0, a real gap develops of 
magnitude U with square root singularities at the band edges. As the temperature increases, the system develops 
substantial subgap DOS which are thermally activated within the ordered phase until T is raised high enough that the 
system enters the normal phase. Plots of the DOS can be found in Ref. |Ht]. Note that the singular behavior occurs for 
one of the "inner" band edges on each sublattice, and that the subgap states develop very rapidly as the temperature 
rises. Furthermore, the DOS on each sublattice is related to the DOS on the other sublattice by a reflection about o = 0. 



INELASTIC LIGHT SCATTERING 



For an electronic system with nearest-neighbor hopping, the interaction with a weak external transverse electromag- 
netic field determined by the vector potential A is described by the Hamiltonian lfl7l,[l8ll : 



he 



2Ti 2 c 2 



(8) 



where the number current operator and stress tensor for itinerant electrons are equal to 

dk„dk 



ja{q) 



^^-dl{k + q/2)d b {k-q/2) and YaJt (q) = £ ^^4 { k + q/2)J b (k-q/2), (9) 



abk 



respectively. Here t a j,(k) are the components of the 2x2 hopping matrix in Eq. (0. The formula for the inelastic light 
scattering cross section derived by Shastry and Shraiman is equal to lfl7l[l8ll 



M ap (q) 



(10) 



It describes the scattering of band electrons by photons with Q. = 0), — (Of and q = kj kf being the transferred 
energy and momentum, respectively, e'^ is the polarization of the initial (final) states of the photons and e«f \ denotes 
the electronic eigenstates. The quantity g{q) = (he 2 /V(Oq) 1 / 2 is called the "scattering strength" with (O q = c\q\, and 
Z = Tr exp(— J3i7) the partition function. The scattering operator M(q) is constructed from both the number current 
operator and the stress tensor; it has both nonresonant and resonant contributions 



M a Hq) i) = (f\Yad9)\i) + L\ 



' (f \j p (k f ) \l)(l\j a (-*,•) | i) (/ \j a (-*,) 1 (/ \jp (k f ) | i) ' 



£/ - £,■ - ft), 



e, - Si + (Of 



(11) 



with the sum I over intermediate states. After substituting into the cross section formula, one obtains three terms in the 
cross section: a nonresonant term; a mixed term; and a pure resonant term (because it is constructed from the square 
of the scattering operator). The components of the cross section can be extracted from the appropriate correlation 
functions (response functions) first calculated on the imaginary Matsubara frequencies and then analytically continued 
onto the real axis Qj^] . Hence, we concentrate on the light-scattering response function ^(^,£2), which is related to 
the cross section and has nonresonant, mixed and resonant contributions: 



R(q,Q.) 



2Kg 2 (k i )g 2 (kf) 

l-exp(-j3£2) 



x(q,&), %(qA) = %N(q&)+%M(qA) + XR(qA) 



(12) 



Inelastic light scattering examines charge excitations of different symmetries by employing polarizers on both the 
incident and scattered light. The Aj g symmetry has the full symmetry of the lattice and is primarily measured by 
taking the initial and final polarizations to be c' = e? = (1, 1, 1, 1, . . .). The Bi 2 symmetry involves crossed polarizers: 
e' = (1,1,1,1,...) and = (—1,1,-1,1,...); while the Z?2g symmetry is rotated by 45 degrees; it requires the 
polarization vectors to satisfy e' = (1,0, 1,0, . . .) and e$ — (0, 1,0, 1, . . .). For Raman scattering (q = 0), it is easy to 
show that for a system with only nearest-neighbor hopping and in the limit of large spatial dimensions, the Ai g sector 
has contributions from nonresonant, mixed and resonant scattering, the Z?i g sector has contributions from nonresonant 
and resonant scattering only, and the Z?2 g sector is purely resonant |[T9il20ll . 

An analysis of the total electronic Raman spectra for the uniform phase of the Falicov-Kimball model has already 
been completed [16]. A full calculation of the nonresonant inelastic light scattering for all q in the CDW chess-board 
phase has also been presented [6]. Here we focus on the total response including the mixed and resonant contributions 
for the Raman scattering (q = 0) in CDW phase. 

MIXED AND RESONANT RESPONSES 

The way we determine the mixed and resonant response functions is as follows: we construct the corresponding 
multi-time correlation function in terms of the generalized polarizations, then perform a Fourier transformation to the 
imaginary Matsubara frequencies and finally analytically continue onto the real frequency axis to extract the response 
function. Such a procedure requires a lot of algebra and analysis tig] ; we do not present all of the details here, but 
instead we summarize the main points. 

The mixed response function is extracted from the multi-time correlation function which is built on three operators: 
one stress tensor and two current operators, as follows 

X ?J M,T 2 ,T3) = (r T 7(Tl)/ (/) (T 2 )J»(T 3 )) . (13) 

The symbol T T is a time ordering operator and the tilde denotes contractions with the polarization vectors [y = 

L e 'a7a ft (<7) e R an d /' — E e ce ja{Tki.f)]- Furthermore, we perform the Fourier transformation from the imaginary 
a/3 ' P a 

times Ti , T2, and T3 to the imaginary Matsubara frequencies iVj,iVf, and zv,- — iVf, As a result, the correlation function 
is represented as a sum over Matsubara frequencies of the generalized polarizations n m+! - , m +i-f m- 

%%f.i(iVi - iVf,iv f , -ivi) = Tj^ [n m _/ >m+/ _/ !m + IT m+1 - >m+i _/ im ] , (14) 

m 

where we introduced the shorthand notation n„,_y .„,+;-/ ,,„ = Tl(ico m — iVf,i(Q m + iVi — iVf,i(O m ) for the dependence on 
the fermionic ico m = inT(2m+l) andbosonic z'V; = ilnTl Matsubara frequencies. In the case of a CDW ordered phase, 
the Feynman diagrams for the generalized polarizations Yl m m+ / are shown in Fig. Q] where we introduce additional 
sublattice indexes a to s. 

Now one has to carefully perform the analytic continuation to the real axis (iVij — > (Qij ± i'0 + , z'V; — iVf — > £2 ± /0 + ) 
and replace the sum over Matsubara frequencies by an integral over the real axis. Then the mixed response function is 
expressed directly in terms of the generalized polarizations 

+00 

= (2^i}2 / d(o\f{(o)-f{(o + a)] 

x Rejn(ct) - (o f + i0 + , co + £2 + i0 + , co - /0 + ) - U(co - a f + z"0+, 03 +Q. - /0 + , co - z0 + ) 

+ n(0) - (Of - z'0 + , co + £2 + ;0 + , co - z0+) - n(<o - co f - z0 + , co + £2 - z'0+ (0 - z0+) (15) 
+ U(co + (Oi + i0 + , co + £2 + z0 + , co - /0 + ) - U(co + co f + z'0 + , co + £2 - z0 + , co - z'0+) 

+ n(a) + cof - /0 + , co + £2 + /0 + , co - z0 + ) - Yl{co + co f -i0 + ,co + Q.- /0+, co - z'0 + ) |, 

where f(co) — 1 /[exp(/3a)) + 1] is the Fermi-Dirac distribution function. 

The mixed contribution is nonzero only for the Ai g symmetry; for other symmetries it vanishes after the summation 
over wave vectors. The next step is to calculate these generalized polarizations. There are two types of diagrams for the 




FIGURE 1. Feynman diagrams for the generalized polarizations of the mixed response function. Due to the properties of the 
dynamic irreducible charge vertex of the Falicov-Kimball model, we have m = in' . 




FIGURE 2. Feynman diagrams for the generalized polarizations of the resonant response function. Both the renormalized and 
bare loop diagrams contribute in all symmetries (A\ g , B\ g and Big)- 

generalized polarizations (see Fig.HJ: the bare loop and the renormalized loop. The reducible charge vertex T ab (shaded 
rectangle in Fig.[T|) is defined from the Bethe-Salpeter-like equation through the irreducible one Y a which is local in 
the DMFT approach on each sublattice [6] and has the same functional form as in the uniform phase lT2lll22ll23ll . Also, 
we used the fact that the total reducible charge vertex is a diagonal function of frequencies for the Falicov-Kimball 
model I2II I22I, l23ll : for other models, where the vertex is no longer diagonal, the analysis is more complicated. The 
final expression for the generalized polarization is too cumbersome to be presented here. 

Similar to the mixed response, the resonant response function is constructed from a multi-time correlation function 
which is built on four current operators 

^/,/,/(T 1 ,T 2 ,T 3) T4) = (7U;%l)/%2)/% 3 )/'V4)) • (16) 

Furthermore, we perform the same formal analytic continuation procedure as for the mixed response function. In 
Fig. 12 we present the Feynman diagrams for the generalized polarizations which contribute to the resonant response 
function. 

For the resonant response function, the analytical continuation onto the real axis is quite complex, but the general 
method remains the same and the final formula for the resonant response function is similar to the mixed response 
function in Eq. ( fT3T > (see Ref. [ 16] for results in the normal phase). In contrast to the nonresonant and mixed response, 
the resonant response contributes to all symmetries. In addition, both the bare and renormalized loops in Fig. [2] 
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FIGURE 3. (a) DOS and (b) different contributions to the Raman responses (co,- = 4) at T = 0.05 and (7 = 2. The solid line 
corresponds to the total response, the dashed-dotted line corresponds to the nonresonant contribution, the dashed line corresponds 
to the resonant contribution, and the dotted line corresponds to the mixed contribution. 
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FIGURE 4. Total Raman response for the A[ g symmetry at T = 0.05 for U = 
Different curves correspond to different incident photon energies ranging from CO, 
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2 on (a) a linear and (b) a logarithmic scale. 
= 1 .5 to co,- = 4.0 in steps of 0.5. 



are present in the resonant response function. In the Z?2 g symmetry we have only the resonant response, in the B\ g 
symmetry we have both nonresonant and resonant responses, and in the A\ g symmetry we have all three responses. 



RESULTS AND CONCLUSIONS 



A detailed analysis of the single particle DOS in the chessboard CDW phase of the Falicov-Kimball model has been 
presented earlier |5J]. Here we show results for Raman scattering in the CDW phase with U = 2 and temperature 
T = 0.05 (which lies below the critical temperature T c w 0.0769). In Fig. [3] we show the DOS (left panel) and the 
different contributions to the Raman response for the A\ g and Z?i g symmetries. In Fig. 0(a) one can see the main 
features of the DOS for chessboard phase at intermediate T: the gap of width U edged by singularities at co = ±tV/2 
and partially filled by subgap states placed at ft) = ±£/2 (£«1 for the case of U ~ 2 and T = 0.05). As a result, the 
scattering spectra displays features (peaks) at frequencies Q. = (U — E)/2, E, (U +E)/2, and U. Such features (peaks) 
were already observed for the optical conductivity ||5J and for the nonresonant Raman and X-ray responses In 
addition, there can also be features at ft), — U, C0j — (U + E)/2, (Oi — E, and ft), ■ — (U — E)/2. 

In the case of B\ g symmetry [Fig. [3] (b)], there is a large peak in the nonresonant response at a frequency equal 
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FIGURE 5. Total Raman response for the B\ g symmetry at T = 0.05 for U = 2 on (a) a linear and (b) a logarithmic scale. 
Different curves correspond to different incident photon energies ranging from ft), = 1.5 to ft),- = 4.0 in steps of 0.5. 
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FIGURE 6. Total Raman response for the Bio symmetry at T = 0.05 for 1/ = 2 on (a) a linear and (b) a logarithmic scale. 
Different curves correspond to different incident photon energies ranging from ft), = 1.5 to ft),- = 4.0 in steps of 0.5. 



to Q. = U that reflects the transitions between the states above and below the gap. The peak has a square root-like 
singularity, that comes from the shape of DOS when there is CDW order and the fact that there is no screening for 
the nonresonant response in this symmetry channel. The smaller peak corresponds to transitions between states of 
the upper (bottom) band and the lower (upper) subgap states at [£l= (U + E)/2]. In addition to the nonresonant 
response, the resonant one has contributions from the renormalized charge excitations (Fig. [2j which reduce the high- 
energy transitions. As was mentioned above, all contributions contribute in the total response for the A\ g symmetry 
[Fig. [3] (b)]. Because of the charge renormalization in the nonresonant, mixed and resonant response functions, the 
total response is smaller than in the B\ g symmetry and there is no square root singularity. 

The total Raman response is presented in set of Figs.[4}|6]for the different symmetry channels and for different en- 
ergies of the incident photons. One can see the main features which were already observed for the optical conductivity 
and nonresonant scattering: four peaks at £2 = (U —E)/2, E, (U+E)/2, and U which correspond to the different in- 
terband transitions. In addition, the resonant and mixed contributions strongly modify the nonresonant response: there 
is a strong enhancement of the scattering when the energy of the incident photon is close to the energy of the interband 
transitions and there is the appearance of additional features (peaks) at the frequencies £2 = (o t ■ — U, O; — (U —E)/2, 
0), — E, and a); — (U + E)/2 as measured from the energy of the incident photon. 

The resonant response is particularly large near the transferred energy £1 w 1 .5 when 0); w 2. By examining results 
on a finer grid of photon frequencies (not shown here) we establish that the resonant profile has a narrow full width at 



half max of much less than 0.1, and a very sharp dependence on o, (the peak height drops by more than three orders of 
magnitude by the time ©, = 2. 1 or 1 .9). In addition, the resonant response in this region does not depend too strongly on 
the symmetry channel of the scattering. Note how similar the curves appear (on a log scale) for the different symmetry 
channels, especially for transferred energies away from U where the gap edge creates sharp features in the B\ % channel. 
Finally, there are joint resonances, as the lower energy peaks do resonate with the large peak, especially for o, w 2. 
Similar resonant effects can be seen for the lower-energy peaks when the incident photon frequency is lower, but they 
are not as dramatic as what happens for the peak near £2 w 1 .5. These resonant effects could be strong signatures of the 
CDW phase in real materials. The resonant enhancements in the normal phase do not produce such enormous peaks 
or have such sharp dependences on the incident photon frequencies; this behavior is arising predominantly from the 
CDW order. 

In conclusion, we have examined the total electronic Raman scattering response for the spinless Falicov-Kimball 
model in the ordered CDW phase. Space limitations allowed us to only consider one value of the interaction and 
temperature, but we see some interesting results, primarily the appearance of a huge resonantly enhanced peak near 
G>i = U that is essentially independent of the symmetry channel. Such a peak could be an important signal for 
experiments on these systems as indicating the appearance of the CDW phase via a direct measurment of the electronic 
charge dynamics. Future work will elaborate on how these features evolve with U and T and will provide additional 
details of the formalism that could not be included here. 
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